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Abstract 

Supersymmetric contributions to time independent asymmetry in J/'4'4' 
process are analyzed in the view of the recent Tevatron experimental measurements. 
We show that the experimental limits of the mass difference AM^^ and the mercury 
EDM significantly constrain the SUSY contribution to B^ — B^ mixing, so that 
sin 2/3s ^ 0.1. We also point out that the one loop SUSY contribution to B'^ J/'4'4' 
decay can be important and can lead to large indirect CP asymmetries which are 
different for different polarization states. These new physics effects in the decay 
amplitude can be consistent with CP measurements in the B^ system. 

1 Introduction 

Recently, CDF and -DO collaborations have announced the observation of CP viola- 
tion in B^ — Bg mixing. The following results, for U^-mixing CP violating phase, 
have been reported [1, 2]: 

2/3, = 0.57 t[J:30 (stat.) 1°:°? (syst.) (DO), (1) 

2/3, G [0.32,2.82] (68%) (CDF). (2) 

These results indicate that the phase Ps deviates more that 3a from the Standard 
Model (SM) prediction [3]. Therefore, the experimental observation of CP violation 
in Bg mixing, along with the Belle and Baber measurement for direct and indirect 
CP asymmetries of B^ decays, open the possibility of probing new physics effect at 
low energy. 



It is a common feature for any physics beyond the SM to possess additional 
sources of CP violation besides the SM phase in quark mixing matrix. In supersym- 
metric extension of the SM, the soft SUSY breaking terms are in general complex 
and can give new contributions to CP violating processes. The SUSY CP violating 
phases can be classified as flavor independent phases, like the phases of the gaugino 
masses and /x term, and flavor-dependent phases, like the phases of the off-diagonal 
^-terms. The flavor independent phases are stringently constrained by the experi- 
mental limits on electric dipole moment (EDAI) of electron and neutron. However, 
the flavor dependent phases are much less constrained. This may imply that SUSY 
CP violation has a flavor off diagonal character just as in the Standard Model. In 
this case the origin of CP violation is closely related to the origin of the flavor 
structures rather than the origin of SUSY breaking [4] . 

The SUSY flavor dependent phases can induce sizeable contributions to direct 
and indirect CP asymmetries of decays [5 7], as in — (f)Ks, Bj^ — rj'Ks and 
Bd Kit which show some discrepancy with the SM expectation. In this paper 
we revisit the supersymmetric contributions to B^ — Bg mixing. We investigate the 
possibility that SUSY may be responsible for the large observed value of Bg mixing 
phase without enhancing the mass difference AMg over the measured value. In 
addition, we analyze the one loop SUSY contribution to B^ J/il^(l> decay, which 
turns out to be important and can lead to a large indirect CP asymmetry. 

The paper is organized as follows. In section 2 we analyze the possible new 
physics contributions to B^—B^ mixing and indirect CP asymmetry of B^ J 
taking into account the constraints imposed by the experimental measurments of 
the mass difference AMg^ and the mercury EDM. In section 3 we discuss the su- 
persymmetric contributions to effective Hamiltonian for AS = 2 and A5 = 1 tran- 
sitions. In section 4 we show that the mercury EDM impose stringent constraints 
on the supersymmetric contribution to the phase > such that the B^ mixing phase 
can not exceed 0.1. In section 5 we analyze the supersymmetric contribution to 
the -Bg J/'4'4' decay. We emphasize that the one loop SUSY contribution to 
i?^ J/ip(p can be important and lead to large indirect CP asymmetries which are 
in general different for different polarization states. Finally, we give our conclusions 
in section 6. 

2 Bg — B*^ mixing and CP asymmetry in B^ J/ip4> 

In the the Bg and Bg system, the flavor eigenstates are given by Bg = (bs) and 
Bg = (bs). The corresponding mass eigenstates are defined as Bl = pBg — qBg 
and Bh = pBg + qBg, where L and H refer to light and heavy mass eigenstates 
respectively. The mixing angles q and p are defined in terms of the transition 
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matrix element Mu = {B^\H^j^='^\B^), where is the effective Hamiltonian 

responsible for AB = 2 transitions: 



P \ M 



12 



(3) 



where we assumed that ATb, < AMb, and ATb, < Tg*"'. The strength of 
mixing is described by the mass deference 



AMb, = Mb„ - Mbj, = 2Re 



= 2|Xi2(B,)|. (4) 



IP 

The decay B^ J/'^(t^ involves vector- vector final states with three polarization 
amplitudes. Therefore, an angular distribution is necessary to separate out the 
three polarizations for a measurement of indirect CP violation without dilution. 
The amplitudes for the decay of B^ f and B^ — > / are given by ^^(/) = 
(/l^i)?-' ISJ) and A\f) = |S;> with 

Here, A, is the polarization index. Therefore, the source of CP violation in decays 
to CP cigcnstates with oscillation arc: oscillation ii q/p ^ 1, decay if p^{f) 7^ 1, 
both oscillation and decay if {q/p,p'^{f)} 7^ 1- The time dependent CP asymmetry 
of -Bg J/'ipfpj for each polarization state A, is given by 

4A _ r>^{mt)^jm-r\B''At)^jm 



CjM cos AMB^t + S),^^ sin AMB.t, (6) 



where Cj^^^ and Sj^^^ represent the direct and the mixing CP asymmetry, respec- 
tively and they are given by 



where r/^ is ± depending on the polarization states. In the SM, the mixing CP 
asymmetry in B^ — J/^^ process is the same for all polarization, to a very good 
approximation, up to a sign. Hence we will omit the polarization index when dis- 
cussing the SM results. We have in the SM, 

sm.2Ps = Sji^^. (8) 

If p{J/%ljct)) = 1, which is the case in SM, then (5s is defined as 2/3^ = arg [M.i2{Bs)\- 
In the SM, the mass difference is given by 

^^If = ^riBmB{BB,Fl;)M^\Vts\''S^{xt). (9) 
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Figure 1: The constraint on R = \Ai\ip/Asm\ in case oi 6 = tt/IO, 7r/4, 7r/2 and 37r/4. 

One may estimate the SM contribution to AMb^ through the ratio AM§^^ / AM§^ , 
where the uncertainties due to short-distance effect cancel: 

AM|f MB,BnJl\V,,r ^''^ 

We can assume that AM§^ = AM^^f ~ 0.507ps-^ Thus, for quark mixing angle 
7 ~ 67 °, one finds AMJ^^ ~ ISps""*^, which is consistent with the recent results 
reported by CDF and DO: 

AMb, = 17.77 ± 0.10(stat.) ±0.07(syst.) (CDF), (11) 

AMb, = 18.53 ±0.93(stat.) ±0.30(syst.) (CDF). (12) 

On the other hand, the SM contribution [p{J /ip(j)) = 1) to the CP asymmetry 
Sj/-^<l, is given by 

Sj,^^ = sin2/?f*^ with /3f^ = arg (^f^) ^ 0{Qm), (13) 

where Vij are the elements of the CKM quark mixing matrix. This result clearly 
conflicts with the experimental measurements reported in Eqs.(l,2). Therefore, a 
confirmation for these measurements would be no doubt signal for new physics be- 
yond the SM. As indicated above, Sj/^^ carries a polarization index corresponding 
to the three final state polarization, however in the SM the mixing induced asym- 
metries are the same( up to a sign) for the three polarizations. 

In a model independent way, the eff'ect of new physics (NP), with p{J /ip(j)) = 1, 
can be described by the dimensionless parameter and a phase 29s defined as 
follows: 

Mf^iBs) '^MfriBsY ^ ' 



Therefore, AMb, = 2\Mf^^ {Bs)\r'^ = AM^^r"^. In this respect, is bounded 
by < AM^^f/AM§f < 1.2. This constrains the ratio between the NP and SM 
amplitudes defined as, R = \Anp/Asm\, as follow: 



1 + i?e*^~^ 



< 1-2 (15) 



Note that for vanishing NP phase, i.e.ONP = 0, one find that R < 0.2. However, for 
Gnp 7^ 0, the constrain on R is relaxed as shown in Fig. 1. It is clear that R can 
be of order one if the NP phase is tuned to be within the range: 7r/2 < 6np < tt, . 

In the presence of NP contribution, the CP asymmetry Bg — ^ J^/V"^ is modified 
and now we have 

Sj/^^ = sin 2/?e// = sin(2/?f ^ + 20,), (16) 

where 

2es = arg (l + Re^""") . (17) 

Therefore, in order to enhance the NP effects, large values of R are required. Now we 
consider the effect of NP that leads to Im[p(J/'0^)] ^ 1. Let us write the amplitude 
as 

A\j/i^^) = A^sMiJm + A^NpiJm, (18) 

and define, 

where 6\ is a weak phase, A is the polarization index, and we have assumed that 
the strong phases in the amplitude ratio cancel. One can now write p(J/V'0) as 

piJ/^l^cP)=e-^<. (20) 

Thus, one obtains, 

-p(J/^(p) = e-2^(^sM+e.+e^)_ (21) 
P 

In this case, the CP asymmetry B^ — ^ is modified and now we have, 

Sj/^^ = sin(2/3f ^ + 2es + 2e\). (22) 

However, as pointed out in Ref.[5], this parametrization is true only when the strong 
phase of the full amphtude is assumed to be the same as the SM amphtude. In fact, 
as discussed in Ref . [8] , the NP strong phase can be different and is generally smaller 
than the SM strong phase thus invalidating the assumption about strong phases 
made in Eq.(19). In general, the SM and NP amplitude can be parameterized as: 

4m = I^Mle^'^'-, ^NP = E I^NPle^'-^e^'^NP, (23) 
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where (5^p arc the strong phases and ^^^P ^^"^ ^IP violating phase. If there is 
one dominant NP amphtude then we can parametrize the NP amphtude as 

^NP = \A^p\e'^»^^e'^i^p . (24) 

Thus, the CP asymmetry Sjj^^ can be approximately written as: 

S)i^^ = sin(2/3f ^ + 20,) + 2r\ cos(2/?f ^ + 20,) sin cos 5^, (25) 

where r\ = I^npMsmI — ^sm ~ '^np- Here A represents the various polar- 

ization states of the vector-vector final state. 

In the SUSY case, considered in this paper, there will be two dominant operators. 
In this case we can write the new physics amplitude as, 

^NP = l^mple^^^'^^e^^f' + \A^iiY'W^^''^ e'^^^^ . (26) 

Now using the result in Ref . [8] , we will neglect the NP strong phase and hence the 
new physics amplitude can be rewritten as an effective single NP amplitude, 

4A _ I 4A \J9^p 
^NP — I^npI^ 



A _ l^iNP I sill ^'^^p + I I sin 6'2Np 

l^lNpl ^INP "I" l^2Npl ^*-*S ^2NP 



tan0- 

I^npI = Vd^mpl sin^iV + \A^^^\ sin^V)' + (l^mpl cos^^^p + I^2NpI cos^V)' (27) 

Hence the expression in Eq.(25) can still be used provided we set the NP strong 
phases to zero. 



3 Supersymmetric contributions to = 2 and = 1 
transitions 

In this section, we analyze the SUSY contribution to the — mixing and 
Bg J/'>P4' decay. As pointed out in Ref. [9], gluino exchanges through AB = 2 
box diagrams give the dominant contribution to B^ — B^ mixing, while the chargino 
exchanges are subdominant and can be neglected. The general H^^^^ induced by 
gluino exchanges can be expressed as 

5 3 
H^/=^ = Y.Ci{^^)Q^i^i) + J2Ci{^i)Q^i^^) + h.c., (28) 

i=l i=l 

where Ci{^), Ci{ij), Qiin) and Qiin) are the Wilson coefficients and operators 
respectively normalized at the scale fj,, with, 

Qi = stj^btsi^X, (29) 
6 



Q2 




(30) 


Qs 


— ^R"L ^R"L-: 


(31) 


Q4 


— ^r"l ^l"r^ 


(32) 


Q5 


— ^r"l ^l"r- 


(33) 



In addition, the operators 0i,2,3 are obtained from <3i,2,3 by exchanging L R. 
The results for the gluino contributions to the above Wilson coefficients at SUSY 
scale, in the frame work of the mass insertion approximation, are give by [10] 



CI 



CI 



CI 



216m| 



24x/6(x) + 66/6(a;) 



216m| 

a? 



23)RLi 



(34) 
(35) 
(36) 



216m| 
'216m| 



2 { [504x/6(x) - 72h{x)] {5i^)LL{5i^)RR - 132/6(x)(4)LR(4)ijL} (37) 



{ [24x/6(x) + 120/6(x)] {5i^)LL{5i^)RR - im^{x){5t^)LR{5i^)RL} (38) 



where x = m-fm- with rrig and rriq are the gluino mass and the average squark 
mass, respectively. The expressions for the functions feix) and f^ix) can be found 
in Ref.[10]. The Wilson coefficients Ci,2,3 are obtained by interchanging the L ^ R 
in the mass insertions appearing in Ci^2.3- 

Note that the mass insertions {^23) L 23) RR may give the dominant contribution 
to the transition matrix element, due to its large coefficient in C|. In order to 
connect Ci{Ms) at the SUSY scale Ms with the corresponding low energy ones 
Ci(/x) with fi ~ 0{mb), one has to solve the RGE for the Wilson coefficients. Also 
the matrix elements of the operators Qi can be found in Ref. [?] . 

Now, we turn to supersymmetric contribution to the amplitude for Bg — ^ J/ip(f). 
It turns out that the gluino exchanges through AB = 1 penguin diagrams gives the 
dominant contributions to this process. The effective Hamiltonian for the AB = 1 
transitions through the penguin process can, in general, be expressed as. 



n 



AB=1 



eff 



^ dOi + CgOg Cidi + CgOg, 



(39) 



i=3 



where 



03 = -sl^^hlcl^^cP^, 

04 = sll^bi^i^^cl, 



(40) 
(41) 
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(42) 
(43) 

(44) 



At the first order in the mass insertion approximation, the gluino contributions to 
the Wilson coefficients Cj^g at the SUSY scale Mg are given by [10] 



Ca{Ms) 
C^{Ms) 
Cq{Ms) 

Cg{Ms) 



^B,{x) + ^B^ix) + ^Pi(x) + lp2(x) 



- ^B2{x) + ^Pi(x) + ^P2{x) 



(45) 



'-Bi{x)-lB2{x)-^Piix) 



rrin 



('^£l)23 ( -M3(x)+3M4(x) ) +(5^^)23^ (^-M3(X)+3M2(X 



1 



The absolute values of the mass insertions (5^^)23, with A,B = {L,R) are 
constrained by the experimental results for the branching ratio of the B Xgj 
decay. These constraints are very weak on the LL and RR mass insertions and the 
only limits we have come from their definition, HS^l Rii)23\ < 1- The LR and RL 
mass insertions are more constrained and, for instance with rrig ~ niq ~ 500 GeV, 
one obtains \{^lr RL)'i3\ ~ 1-6 x 10~^ [7, 10]. Note that, although, the LR{RL) mass 
insertion are constrained severely their effects to the decay are enhanced by a large 
factor my/niiy as can be seen from the above expression for Cy{Ms)- 

In this respect, the phase of {^'[^1)23, (^ll)23 ^^'^ ('^r_r)23 t^is relevant CP 
violating phases for our process. In the next section, we discuss possible constraints 
imposed on these phases by the mercury EDM. 



4 Mercury EDM versus large Bg — Bg mixing phase 

It has been pointed out [12, 13] that large values of (^23)^^? enhance the chromo- 
electric dipole moment of the strange quark which is constrained by the experimental 
bound on the EDM of mercury atom Hg. In this section we show that the Hg EDM 
imposes a constraint on l^n^ii^LL) 23 {^rr) 23], which may limit the supersymmetric 
contribution to the B'^^ — B^ mixing. 

In the chiral lagrangian approach, the mercury EDM is given by [12] 

dHg = -e (d^ -d^- 0.012df ) x 3.2 x 10"^ (46) 
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The chromelectric EDM of the strange quark df^ is given by 



jC _ 9s 



m 



(47) 



where x = m?/m|. For = 500 GeV and x = 1, the experimental limit on Hg 
EDM leads to the following constraint on {623) lr- 

Iia{S^2) LR< 5.6 X 10'^ . (48) 

The mass insertion ((^22)^^ may be generated effectively through three mass inser- 
tions as follows: 



(5; 



22) LR 



i^23)LL{S33)LR{^32)RR^ 



(49) 



where {6i^)LR ^ ^ C'(10~2)_ Therefore, the Hg EDM imposes the 

d 

following constraint on the LL and RR mixing between the second and the third 
generations: 



Im 



< 5.6 X 10" 



(50) 



(4)ll(4)Lj 

If one assumes that (^23)ii ~ with negligible weak phase, then he gets the 
following bound on the {S2s)rr mass insertion: 



\{^23)rr\ sin 



(51) 



Therefore, in case |((523)^fil ~ 0(0.01), the associated weak phase is essentially 
unconstrained. However, if |(523)-R-r| ~ 0(0. 1), the the weak phase is constrained to 
be of order 0.1. In both cases, this will limit the SUSY contributions to the Eg — Bg 
mixing phase. 

We start our analysis for SUSY contribution to sin 2/3s by assuming that Bg — Bg 
mixing may receive a significant SUSY contribution, while the decay of Bg — J/i/jcf) 
is dominated by the SM. Therefore, we have lm[p{J/ilj(j))] = and the induced CP 
asymmetry is given by Sj/^^ = sin(2/3f^ + 26 g). As an example for the SUSY 
contribution, we consider nig = 500 GeV and x = 1, which leads to the following 
expression for R = \Mf^^^ /Mf^^ [9]: 



R 



1.44 



('^23)iL + ('^23)rR 



■d \2 



+ 27.57 



+ (4) 



RL 



44.76 



{Sis)LR{S23)RL 



- 175.79 



(^23)ii('^23)fii? 



(52) 



Prom this equation, it is noticeable that the dominant contribution to the Bg — Bg 
mixing is due to the mass insertions {^23) l 1(^23) RR- 

If one assumes that (^23)^^ induced by the running from the high scale, where 
left-handed squark masses are universal, down to the electrweak scale, then one 
finds (^23)^'^' ~ ~ 0.04. With a small source of non-universality in the right- 
handed squark sector, one can easily get ((^23)^^ of order 0(0.1). Therefore, one 
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Figure 2: The — mixing phase as function of the arg[(523 radians) for 
1(5^3)^^1 = 0.025,0.05 and 0.1. 

gets i? ~ 0.7. However in this case, the Hg EDM implies that: arg[(5^j^)23] < 0.1, 
which hmits significantly the SUSY effect for enhancing sin2/5s. 

In Fig. 3, we present our results for the -B^ — B'^ mixing phase 2/3^ as a function 
of eLTg[{S2s) rr] for \{S23)rr\ = 0.025,0.05 and 0.1. At these values the ratio R is of 
order < 0.17, 0.35 and 0.7 respectively. As can be seen from this figure, the values of 
Bg mixing phase, which are consistent with the Hg EDM constraints, are typically 
of order < 0.1. Therefore, one concludes that the SUSY contribution through the 
Bg — Bg mixing implies limited enhancement for sin2/3s and thus cannot account 
for the new experimental results reported in Eq.(l,2). Moreover, a salient feature 
of this scenario with large RR mixing is that it predicts a reachable mercury EDM 
in the future experiments. 

5 SUSY contribution to J/ip4' decay 

In this section we will consider SUSY contribution to the decay B'^ J/tpcl). How- 
ever, let us discuss the complexities in analyzing new physics effects in the decay 
amplitude for vector- vector final state[14]. 

Consider aB V1V2 decay which is dominated by a single weak decay amplitude 
within the SM. This holds for processes which are described by the quark-level 
decays b — > ccs which is the underlying quark transition in B^ — > J/tjjcj). In this case, 
the weak phase of the SM amplitude is zero in the standard parametrization [15]. 
Suppose now that there is a single dominant new physics amplitude, with a different 
weak phase, that contributes to the decay. This indeed will be the case for the SUSY 
contribution to —>■ J/Tp(j). The decay amplitude for each of the three possible 
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helicity states may be written as 

= Amp(5 ^ ViV2)x = OAC^^A + b^e^'^e'^x , 
Ax = Amp{B ^ ViV2)x = axe'^l + bxe-'^e< , (53) 

where ax and bx represent the SM and NP amplitudes, respectively, (/) is the new- 
physics weak phase, the 5^'^ are the strong phases, and the helicity index A takes the 
values {0, II, ±}. Using CPT invariance, the full decay amplitudes can be written as 

A = Amp{B V1V2) = Aogo + A\\g\\ + i A±g^_ , 

A = Amp(B ^ V1V2) = Aogo + - i A±g± , (54) 

where the gx are the coefficients of the helicity amplitudes written in the linear 
polarization basis. The gx depend only on the angles describing the kinematics [16]. 
Eqs. (53) and (54) above enable us to write the time-dependent decay rates as 

r{B^,{t) ^ V1V2) = e-r* ^ (Aa, ± J^xa cos(AMt) ^ pxa sm{AMt))gxga • (55) 

A<fT 

Thus, by performing a time-dependent angular analysis of the decay B(t) — > V1V2, 
one can measure 18 observables. These are: 



Aaa 


= l{\Ax\' + \Ax\'), 


Saa 


= \(.\A,\'-\Axn 


A±i = 


= -lm{A^A*-A^A*), 


ho-- 


= Re{A\\A*+A\\A*). 


S±i = 


= -lm{A^A*+A±A*), 


S||o ■ 


= Re(^l|A*-A||i*) 


P±i = 


R^^[AlAi+A*A^]), 


P±± 


= lm(^^lA^), 


P\\o = - 


-lm{^[AlAo+A*A\\]), 




-lm(UlA,), 



where i = {0, ||}. In the above, q/p is the weak phase factor associated with B^-Bg 
mixing. For Bg meson, q/p = exp(— 2z/?s). Note that Ps may include NP effects 
in Bg-Bg mixing. Note also that the signs of the various pxx terms depend on the 
CP-parity of the various helicity states. We have chosen the sign of poo and p\\\\ to 
be —1, which corresponds to the final state J/ip<p. 

Not all of the 18 observables are independent. There are a total of six ampli- 
tudes describing B V1V2 and B V1V2 decays [Eq. (53)]. Thus, at best one 
can measure the magnitudes and relative phases of these six amplitudes, giving 11 
independent measurements. 

The 18 observables given above can be written in terms of 13 theoretical pa- 
rameters: three ax^s, three bx^s, (5s-, (p, and five strong phase differences defined by 
dx = S\ — 5"^, Ai = 5'j_ — 5f . The explicit expressions for the observables can be 
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Figure 3: SM tree level (left) and SUSY one loop (right) contributions to — > J/'ip4> 
decay. 

found in Ref [14]. In the presence of new physics, one cannot extract the phase 
Ps- There are 11 independent observables, but 13 theoretical parameters. Since the 
number of measurements is fewer than the number of parameters, one cannot ex- 
press any of the theoretical unknowns purely in terms of observables. In particular, 
it is impossible to extract Ps cleanly. 

In the absence of NP, the bx are zero in Eq. (53). The number of parameters is 
then reduced from 13 to 6: three aa's, two strong phase differences (Aj), and Ps- It is 
straightforward to show that all six parameters can be determined cleanly in terms 
of the observables. This is exactly what is done in the experimental measurements 
to measure (3s, the value of which appears to be inconsistent with the SM. This 
might indicate new non SM phase in Bs mixing or NP in the decay amplitude in 
which case the general angular analysis in Eq. (55) should be used. In the presence 
of NP, the indirect CP asymmetries for the various polarization states are not longer 
the same as it is in the SM( up to a sign). 

In this section we will consider the scenario where SUSY gives significant con- 
tribution to both B° - B° mixing and the decay of Bg — > J/ip(f). In this case, the 
induced CP asymmetry is given by Eq.(25). As shown in Fig. (5), the SM the decay 
of B^ — > J/i^xj) takes place at tree level through the b ^ c transition. While the 
dominant SUSY contribution to this decay is given by the one loop level of gluino 
exchange for b ^ s transition. It is interesting to note that the SM amplitude is 
proportional to Gp x V^cKs ~ 10~^, while the SUSY amplitude is given in terms of 



it can be important relative to the SM one. In this respect, it is important to con- 
sider the impact of this contribution on the induced CP asymmetry Sj^^^, as the 
phase of the mass insertion ((5^^)23 is not constrained by EDM. 

Let us now write down the SM and SUSY contribution to Bg{p) — J/ijj{ki,ei)<p{k2, €2), 
where we have labelled the momentum and polarization of the final state particles. 
To proceed with our calculation, we will first specify the momentum and polariza- 




Therefore, although SUSY contribution is a loop level. 
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tion vectors of the final-state particles. We will work in the rest frame of the 
meson. We define the momentum and polarization of the vector (p meson as [17] 



4(0) = —{-k,0,0,E^) 

TTlfj) 

= -^(0,Tl,-i,0) , (57) 
The momentum and polarization vectors for J/^|J are defined as, 

= ^(0,Tl,-^,0) , (58) 

The general amplitude for -B^(p) — > J/ip{ki,ei)(p{k2,£2), can be expressed as[18], 

A = ael-e*2 + -^{p ■ e*^){p ■ e^) + i^e^,p^p^'q''el"er , (59) 
rrip, mi^ 

where q = ki — k2- For angular analysis it is useful to use the linear polarization 
basis. In this basis, one decomposes the decay amplitude into components in which 
the polarizations of the final-state vector mesons are either longitudinal (^o)) or 
transverse to their directions of motion and parallel (Ay) or perpendicular {A±) to 
one another. One writes [19, 20] 

- - 1 - i - 

A A ^ A -*T -*T A -*T ^, -*T t ar\\ 

A = AQei -£2 - -j=A\\ei -eg ~ -P, (60) 

where p is the unit vector along the direction of motion of V2 in the rest frame of 
Vi, = -p, and ef'^ = ^1 — ^t^P- Aq, , A± are related to a, b and c of Eq. (59) 
via 

A.=V2a, Ao = -ax-'^bix'-l) , = 2^2 ^^cv^^^ , (61) 

where x = ki ■ k2/{mim2)- (A popular alternative basis is to express the decay 
amplitude in terms of hclicity amplitudes A\, where A = 1,0,-1 [19,21]. The 
helicity amplitudes can be written in terms of the linear polarization amplitudes via 

A±i = (A|| lb A_|_)/-v/2, with ^0 the same in both bases.) 

We will now proceed to calculate the polarization dependent CP asymmetry 
given in Eq. 25. We will use factorization to calculate the ratio r\ = |A^p/Agp^|. 
In factorization there are no strong phases and we will keep them as a free unknown 
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parameter in the expression for Sj^_^^ in Eq. (25). The ampUtude for the process 
Bsip) J/i'{ki,£i)4>{k2,<^2), in the SM, is given by, 

A[Bs ^ Jm = -^XLj,^, (62) 



with 



X = VcbV*ga2- VqbV*s{al + al + a'^j + a? 



9J ' 
q=u,c,t 

Lj/^ = mj/^gj/^el^ {(l)\sjfj,{l-j5)b\Bs) , . (63) 

where a2 = C2 + -^ and for i > 2, ai = Q + t^, with Cj being the Wilson's coefficient. 
Here gj/^ is the J/ip decay constant defined in the usual manner. 

We can simplify X using several facts. First 02 is much larger than a* with 
i = 3, 5, 7, 9 [22]. Second, in the penguin contributions in Eq. (63) we have included 
the rescattering contributions from the tree operators. However these are small and 
the contributions a^'^ and a^'^ due to perturbative QCD rescattering vanish because 
of the following relations, 

4$ = -ct;,/N, = Pr/Nc, (64) 

where is the number of color. The leading contributions to P* are given by: 
PI = (^)ci(^ + G{mi, ix,(f')) with i = u,c. The function G(m, /U,g^) is given by 

G{m, n, q^) =4 r x{l - a:)ln '"' ~ ""^l ~ "^^^^ dx . (65) 
Jo IJ- 

The rescattering via electroweak interactions are given by [23] 

Cr:9 = Pe^', cZo = (66) 
with Pg = {^^){NcC2 + ci)(^ + G{mi, fi, q^))- These contributions are again much 
smaller than the dominant tree contributions and can be neglected. 

In light of the above facts we can conclude that the dominant contribTitions in 
X in Eq. (63) come the tree level term where ci = 1.081 and C2 = —0.190 are the 
relevant Wilson coefficients [22]. This leads to 

X « ye6K*s«-2 = 0.17^6 (67) 

The matrix elements in Eq. (63) above can be expressed in terms of form factors. 
This can be done as follows [24]: 



{V2{k2)\ct^^,b\Bs{p)) = i- ^^,e^upaP''k'2^7 

[iriB + m2) 

{V2{k2)\q'l^,j5b\B{p)) = {mB + m2)Af\r^) 



n2\ 



^2-r 



TUB + m2 



{Pij. + k2^) 



rrin — mr, 



+ 2im2'-fr^A^o\r^) > (68) 
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where r = p — k2, and y(^), A^^\ A^^ and A^'^ are form factors. 
Using Eq. (68) in Eq. (63) one obtains, 



asM 



-^■mj/^i,gj/^i,x{mB, + m^)Af\rr?j/^)X 

bsM = 

Gf tub, 
csM = — T^mj/^gj/^- ■ — 



-mB^V('\m]/^)X. 
<t>) 



(69) 



Let us turn now to the SUSY contribution. We will consider only the dominanat 
chromomagnetic operators. The gluon in these operators can split into a charm 
quark pair, thereby contributing to 6 — > sec. We begin with a discussion on the 
matrix elements of the chromomagnetic operators Og and Og. These are given by. 



{J/i;4>\Oa\Bs) 



<Og> 

asUib 



{J/i;cf>\Og\Bs) = <Og> 



2 



\Bs) 
(70) 



where is the momentum carried by the gluon in the penguin diagram. In our 
case q'^ coincides with the four momentum of the J/V'- 
After a color fierz we can write the operator Og as, 

(5a7/.^(l + 75)6a) {m^s^) + .. 





= Y9 


\--{ 




= Y9 


I AT J 











In the above we have only retained terms that contribute to the decay Bs{p) 
J/ip{ki, ei)4>{k2, S2) In factorization, after using equation of motion, we can write 
the matrix element of Og as. 



Ti 



OgYg 



2 . 



L 



■mj/^gj/^e*/ {(t>\ sj^{l - 75)6 \Bs) , 



T2 = CgYg — 



2 r. ■ 
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R 



mj/^gj/^el^ {4>\ s-f^il + 75)6 \Bs) , 

2£l ■ k2 



nib 



2 ■ 



mj/^gj/^ (01 s(l + 75)& \Bs) ■ 



(71) 



In the above m^^t are the strange and the bottom quark masses. 

In the above equation it is clear that T2 is suppressed relative to Ti by ^ and we 
will neglect it. Prom the structure of the polarization vectors in Eq. (57), it is also 
clear that the it polizations do not contribute to T^. Hence for the it polarizations 
we can obtain a clear prediction for r\ defined below Eq. (25), as the form factors 
and other hadronic quantities cancel in the ratio. 

For the matrix element of the operator Og, focussing only on the transverse 
amplitues we can write, 



<Og> 



2 „ 



Rj/^ = ■mj/^gj/^e*/ {(f>\s'yij,{l + 'y5)b\Bs) , 



(72) 



Hence again focussing only on the transverse amplitudes we can write, using 
Eq. (68) in Eq. (71) and Eq. (72), 



^susy 



^susy 



Y = 



Y 



-^rnj/^gj/^{mB, + m^)Af\m^j/^){Y - Y) 
Gf 

--^rnj/^gjl^ 
\/2C, 



{tub, + rriA 



mBj^''\m''j,^){Y + Y) 



Gf_ 



Ya 



Yn 



(73) 



Combining the SM and SUSY contributions we can now compute. 



''a — I^npMsmI — 



SMi 



X ^ 
{Y + Y) 
X ' 



(74) 



Using the values of Vcb and Vcs from Ref [15] we obtain X k. 0.0069. Futhermore 
with m~ = m~ = 500 GeV, mb{inb) = 4.5 GeV, we obtain, 



Y 
Y 



= 2.1315(<^iK)23 
= 2.1315(4^)23 



= 0.0477(5^^)23 



(75) 
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We can then write, using Eq. 74, 



;(lfe)23|)^ + ilKMr - 2|fe)23||fe)23| COS{0LR - 0Rl) 

0.01 



rA ~ 0.07 — (,76) 

where 9lr and 9rl are the phases of ((^fij)23 and ((^^i)23. We will set |((^fij)23| = 
I ((5^^)23 1 = 0.01 and we can then now consider the following cases: 



Case a {S'lji)23 = {^rl)23- In this case we obtain. 



sLi,. = sin(2/3f ^ + 20, 



Sj/^^ = sin(2/3f^ + 20,) + O.28cos(2/3f^ + 26',)sine^pCOS(5^. (77) 

If we neglect the contribntion from mixing then Sj^^^ can reach a value of upto 0.3 
for sin^jijp ^ 1 and cos (5-*- ~ 1. 

Case b {S'lji)23 = —{Srl)23- In this case we obtain, 

-^1/^0 = sin(2/?f ^ + 29s) + 0.28 cos(2/?f ^ + 29s) sin ^Jp cos (5" 

Sj/^^ = sin(2/3f^ + 20,). (78) 

Again, if we neglect the contribution from mixing then Sj^^^ can reach a value of 
upto 0.3 for sin^jljp 1 and cosj" ~ 1. Finally, we can consider the case where 
either {S'Ir)23 or (<^^^)23 is zero. For the case ((^f ^2)23 7^ 0, (5^^)23 = we obtain, 

S]/^^ = sin(2/3f ^ + 29s) + 0.14 cos(2/3f ^ + 29s) sin 4p cos 5" 
For the case {S'Ir)23 = 0, {Srl)23 / we obtain. 



.Sj/^^ = sin(2/3f^ + 20,)+O.14cos(2/3f^ + 20,)sin4pCOs5^ (79) 



= sin(2/3f ^ + 29s)- 0.14cos(2/3f ^ + 20,) sin 0|p cos (j" 

= sin(2/3f^ + 20,)+O.14cos(2/?f^ + 20,)sin0^pCOS(5^ (80) 



Now one may wonder how NP in 6 ^ sec transitions affect CP measurements 
in the system. Let us first consider the indirect CP asymmetry in the golden 
mode Bd J/ipKs- Note this is a vector-pseudoscalar decay and so the strong 
phases involved here can be quite different from the ones involved in vector-vector 
decays. In other words NP effects in different final states can be very different. More 
interestingly, it can be easily checked that for case b in Eq. (78) the contribution 
to the indirect asymmetry in B^ J/ipKs cancels. However, the indirect CP 
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asymmetry in the vector-vector mode does not cancel for all polarization states. In 
other words the range of NP effects obtained in the decay Bg J/ip(p are consistent 
with sin 2/3 measurements in — J/ipKg [25-27] for the various reasons discussed 
above. 

The decay B^ — J/il^K* is related to B^ — J/ilj<p by SU{3) flavor symmetry. 
Hence we should potentially see NP effects in B^, — J/ijjK* , up to SU{3) breaking 
effects. The CP measurements in this decay are not yet precise [25] and hence this 
decay also is an ideal place to look for new physics effects in the decay amplitude. 

5.1 Summary 

In summary, we have analyzed the SUSY contribution to B^ — B^ mixing in light 
of recent experimental measurement of the mixing phase. We showed that the 
experimental limits of the mass difference AMb^ and the mercury EDM constrain 
significantly the SUSY contribution to B^ — Bg mixing, so that sin2/3s < 0.1. We 
then studied the the one loop SUSY contribution to B^ J/4'4> decay and found 
that new physics contribution to the decay amplitude can lead to significant indirect 
CP asymmetries which are in general different for different polarization states. 
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